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Abstract 

Mutations in the F-box only protein 7 gene (FBX07), the substrate-specifying subunit of SCF E3 ubiquitin ligase complex, 
cause Parkinson's disease (PD)-15 (PARK15). To identify new variants, we sequenced FBX07 cDNA in 80 Taiwanese early 
onset PD patients (age at onset <50) and only two known variants, Y52C (c.155A>G) and Ml 151 (c.345G>A), were found. 
To assess the association of Y52C and Ml 151 with the risk of PD, we conducted a case-control study in a cohort of PD and 
ethnically matched controls. There was a nominal difference in the Y52C G allele frequency between PD and controls 
(p = 0.045). After combining data from China [1], significant difference in the Y52C G allele frequency between PD and 
controls (p = 0.012) and significant association of G allele with decreased PD risk (p = 0.017) can be demonstrated. Upon 
expressing EGFP-tagged Cys52 FBX07 in cells, a significantly reduced rate of FBX07 protein decay was observed when 
compared with cells expressing Tyr52 FBX07. In silico modeling of Cys52 exhibited a more stable feature than Tyr52. In cells 
expressing Cys52 FBX07, the level of TNF receptor-associated factor 2 (TRAF2) was significantly reduced. Moreover, Cys52 
FBX07 showed stronger interaction with TRAF2 and promoted TRAF2 ubiquitination, which may be responsible for the 
reduced TRAF2 expression in Cys52 cells. After induced differentiation, SH-SY5Y cells expressing Cys52 FBX07 displayed 
increased neuronal outgrowth. We therefore hypothesize that Cys52 variant of FBX07 may contribute to reduced PD 
susceptibility in Chinese. 
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Introduction 

Parkinson's disease (PD) is attributable to environmental factors, 
whereas evidence has shown that genetic factors are also playing 
an important role. Molecular genetic studies have identified 
thirteen genes linked to rare dominant or recessive monogenic 
forms of PD: SNCA, Parkin, PINK1, DJ-1, LRRK2, ATP13A2, 
VPS35, PLA2G6, FBX07, EIF4G1, SYXJ1, DNAJC6, and 
DNAJC13 [2,3]. 

Recently, the F-box only protein 7 (FBX07) mutations have 
been identified in several families with early-onset parkinsonism 
and pyramidal tract signs. Homozygous R378G missense muta- 
tion in an Iranian kindred, homozygous nonsense mutation 
(R498X) in an Italian family, a Pakistan family and a Turkey 
family, and compound heterozygous mutations (IVS7+1G/T and 
T22M) in a Dutch family are unambiguously responsible for the 
autosomal recessive, early-onset, parkinsonian-pyramidal syn- 
drome [4,5,6]. The presynaptic nature of the parkinsonism in 
the families reported by Di Fonzo et al. is shown by the dramatic 
abnormality of DaTSCAN SPECT, the beneficial effect of 
levodopa, and the presence of levodopa-induced dyskinesias, 
suggesting that FBX07 mutations may be the potential genetic 



causes of early-onset Parkinson's disease (EOPD) or familial PD 
[5]. Two missense substitutions, p.Ile87Thr and p.Asp328Arg, in a 
single heterozygous state, were found in two EOPD patients in 
Taiwan [7]. Although no pathogenetic mutations in the FBX07 
gene were detected in 135 Chinese early-onset parkinsonism 
patients, the PD patients showed a trend toward decrease in Y52C 
G allele frequency compared with the controls [1]. 

FBX07, a member of the F-box-containing protein (FBP) 
family, encodes a protein of 522 amino acids consisting of several 
discrete domains: the ubiquitin-like domain, cyclin-dependent 
protein kinase 6 (CDK6) binding site, FBX07/PI31 domain, F- 
box motif, proline-rich region, and R(ar)DP motif [5,8] . Through 
the interaction between the F-box and the Skpl protein, FBPs 
become part of SCF (Skpl-Cullinl -F-box protein) ubiquitin ligase 
complexes, and play roles in ubiquitin-mediated proteasomal 
degradation (review in [9]). F-box proteins recruit a large number 
of diverse substrates to SCF complexes and allow for their 
ubiquitination [10]. FBX07 promotes ubiquitin conjugation to 
TRAF2 (a member of the tumor necrosis factor receptor 
associated factor protein family with ubiquitin ligase activity) and 
cIAPl (an apoptosis inhibitor possessing ubiquitin ligase activity), 
resulting in decreased receptor-interacting protein 1 (RIP1) 
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ubiquitination and lowered NF-kB signaling activity [11,12]. 
Strong evidence has shown that NF-kB induced neuroinflamma- 
tion may be involved in development of PD [13,14]. Therefore, 
FBX07 may play a role in protecting neurons from PD process. In 
contrast, the mutations or variations that change the function, 
expression, or stability of FBX07 may confer PD risk to the 
subjects. 

Results 

Mutation/variant analysis of FBX07 

Since most of the mutations found in FBX07 result in truncated 
FBX07 protein, amino acid replacement or multiple aberrant 
frame-shift splice variants, we sequenced FBX07 cDNA (Table 1) 
instead of genomic DNA in a cohort of ethnic Chinese patients 
with EOPD in Taiwan to identify previously undetected variants 
in EOPD cases of Chinese origin, followed by a case-control study 
for the identified variants. The cDNA samples contained only 
DNA sequences from genes that were transcribed into RNA. Thus 
FBX07 cDNA fragments from 80 EOPD patients were amplified 
for sequence analysis. However, only two known substitutions that 
caused changes in the peptide sequence were identified: a 
c. 155A>G substitution leading to an amino acid change from 
tyrosine to cysteine in position 52 (Y52C) in one EOPD patient 
(heterozygote) and a c.345G>A substitution resulting in a 
methionine to isoleucine change at amino acid position 1 1 5 
(Ml 151) (rsl 1 107) in 74 EOPD patients (40 homozygotes and 34 
heterozygotes) (amino acid number according to NM_012179) 
(Fig. 1A). The two reported [1] variants were confirmed using 
PCR-restriction fragment length polymorphism (RFLP) method 
(Fig. IB). Both Y52C and Ml 151 are not evolutionarily conserved 
in the known mammalian homologues of the FBX07 protein 
(Fig. 1C). 

Case-control study of Y52C and Ml 1 51 

A case-control study in a cohort of PD patients (n = 5 1 6, 80 
EOPD patients included) and ethnically matched controls 
(n = 516) was conducted to assess the association of Y52C and 
Ml 151 with the risk of PD (Table 2). All genotype frequencies 
were confirmed to be in the Hardy-Weinberg equilibrium. There 
was no statistically significant difference (/>>0.025) in genotype or 
allele distribution between patients and controls for both single 
nucleotide polymorphisms (SNPs) examined, after correction of 
multiple SNP testing. However, for Y52C, the frequency of AG 
genotype (0.8% vs. 2.3%, p = 0.044) or G allele (0.4% vs. 1.2%, /> 
= 0.045) was notably lower in PD patients than the controls. Y52C 
AG genotype or G allele demonstrated a trend toward decrease in 
risk of developing PD (odds ratio: 0.33, 95% confidence interval: 
0.09-0.95, p — 0.055-0.056). Analysis combining our patient and 
control subjects as well as the population in Luo's study [1] yielded 
results of statistically significant difference in genotype (0.9% vs. 
2.8%, p = 0.012) and allele (0.5% vs. 1.4%, p = 0.012) distribution 
between patients and controls. The negative association of the 
Y52C AG genotype or G allele with PD was significant (odds ratio: 
0.32-0.33, 95% confidence interval: 0.12-0.77, /> = 0.0 16-0.0 17). 
The identified one EOPD and three late-onset PD patients 
carrying Y52C presented with asymmetrical tremor, rigidity, and 
bradykinesia without pyramidal signs, all of whom had a good 
response to anti-parkinsonian medication. 

FBX07 expression analysis 

Since the case-control study suggests that FBX07 Y52C G allele 
might be a potential protective factor, we cloned the polymorphic 



FBX07 cDNA, which was then expressed in HEK-293T cells to 
investigate the functional consequences. 

The common cellular abnormality found in the PARK 15 
patients from the Dutch and Italian families is the depletion of the 
FBX07 isoform 1 (NM_012179), which normally is located in the 
cell nucleus [15]. Forty-eight hours after transfection of EGFP 
tagged FBX07 isoform 1 constructs, cells were analyzed by 
fluorescent microscopy. Although Cys52 FBX07 protein dis- 
played nuclear and cytosolic staining pattern similar to Tyr52, a 
significantly stronger green (FBX07 fusion protein) relative to blue 
(nuclei staining) fluorescence signal was observed in Cys52 
FBX07 cells (1.60 vs. 2.57, /> = 0.015; Fig. 2A). To further 
examine the transiently expressed FBX07-EGFP fusion protein, 
protein blotted with FBX07 antibody was performed. As shown in 
Fig. 2B, FBX07-EGFP fusion protein in the expected size range 
for Tyr52 and Cys52 constructs was observed. However, the 
protein expression level of Cys52 FBX07 was increased compared 
with the Tyr52 FBX07 (21 1 %,p = 0.016). The stability of Cys52 
variant was further examined by a cycloheximide (200 ug/ml) 
chase experiment. While the Tyr52 protein was degraded to 68%, 
18%, 11%, 9% and 7% left after 6, 12, 24, 36, and 48 hr of 
protein synthesis blocking, reduced rates of decay were observed 
for Cys52 variant (90%, 78%, 72%, 52%, and 21% remained, 
respectively) (Fig. 2C). 

Homology modeling of Cys52 FBX07 

To understand the structure-based information of Cys52 variant 
in FBX07, homology modeling of Tyr52 and Cys52 FBX07 was 
performed. After energy minimization, the modeled structures for 
Tyr52 (WT) and Cys52 (Y52C) were shown in Fig. 3. The 
potential energy of Tyr52 and Cys52 mutant was —2463.854 and 
-2471.736 kcal/mol, indicating Cys52 FBX07 exhibited a more 
stable feature than Tyr52 FBX07. According to hydrogen-bond 
(H-bond) computing analysis, the H-bond interaction of Tyr and 
Cys52 was shown. In the Tyr52 model, Tyr52 did not form any 
H-bonds with adjacent residue. On the other hand, H-bond 
formed by the Cys52 with the Asp54 causes decrease of the local 
energy. 

FBX07 regulating TRAF2 abundance through 
ubiquitinating TRAF2 

Through binding and mediating ubiquitin conjugation to 
TRAF2 and c.IAPl, FBX07 was identified as a negative regulator 
of NF-kB signaling [12]. To assess Cys52's effect on TRAF2 
abundance, Tyr52 or Cys52 FBX07 cDNA plasmid was 
transfected into HEK-293T cells and protein was blotted with 
TRAF2 and FBX07 antibodies (Fig. 4). Compared with the 
endogenous FBX07 level, Tyr52 or Cys52 cDNA transfection 
significantly increased FBX07 abundance (3.11-7.11 folds, p< 
0.001). The FBX07 level of Cys52 was significantly higher than 
that ofTyr52 (6.15-7.11 vs. 3.1 1-4.01, p = 0.001). Accompanying 
that, TRAF2 protein was significantly decreased in Tyr52 or 
Cys52 FBX07 cells (0.69-0. 0.89, p = 0.008-0.001). In Cys52 
cells, the TRAF2 expression level was significantly lower than that 
of Tyr52 cells (0.69-0.70 vs. 0.80-0.89, p = 0.010). 

Furthermore, to examine whether Cys52 exerts effects on 
TRAF2 ubiquitination, proteasome inhibitor MG-132 (5 uM) was 
applied to Tyr52 or Cys52 FBX07-EGFP expressing HEK-293T 
cells and co-immunoprecipitation was performed (Fig. 5A). The 
results showed a significant increase in FBX07-TRAF2 interac- 
tion (367% vs. 100%, p= 0.019, Fig. 5B) in Cys52-expressing cells 
as compared to Tyr52-expressing cells. Moreover, a significant 
increase in TRAF2 ubiquitination was observed in Cys52 cells as 
compared to Tyr52 cells (132% vs. 100%, p = 0.028) and in Tyr52 
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Y52C (TAC/TGC) 
A A C TA/G C A A G 




M115I (ATG>ATA) 
A G CA TG/ACA G 




B 



Pst\ RFLP (CTGCAG) 
AA AA AA AA AG AA 




Pael RFLP (GCATGC) 
AA AA GG GG AG AG 




Human 
Pan troglodytes 
Bos taurus 
Mus musculus 
Rattus norvegicus 



52 

NTRFTITLNYKDPLTGDEE 
NTRFTITLNYKDPLTGDEE 
DARFAITLNNKDALTGDEE 
DTRFAITLNNKDALTGDEE 
DTRFAITLNNKDALTGDEE 



115 

LATSSNQTSMQDEQPSDSF 
LATSSNQTSMQDEQPSDSF 
LATSSSQSNIQDAQLHDSL 
LAATPSQTNIPDEQGTDSS 
LAATSSQANIPDEQGSDSS 



Figure 1 . Variant identification and amino acid sequence alignment. (A) Chromatograms of direct cDNA sequencing of Y52C and Ml 1 51. (B) 

Restriction analysis of Y52C and Ml 151. On agarose gel, Y52C creates new Psfl restriction site by mismatch PCR and leads to an additional 222 bp 
band, whereas M1 151 loses restriction by Pael and leads to 479 bp band. (C) Evolutionary conservation of the regions of FBX07 Y52C and M1 151 using 
the program Vector NTI. 
doi:1 0.1 371 /journal.pone.01 01 392.g001 



cells as compared to vector-transfected cells (100% vs. 80%, 
/i = 0.0 1 5) (Fig. 5B). These results suggest that FBX07 may 
decrease TRAF2 abundance through promoting TRAF2 ubiqui- 
tination. 

Suppression or over-expression of FBX07 affecting cell 
survival upon MPP + treatment 

To gain further insight into the possible role of FBX07 protein 
involved in PD, the mitochondrial inhibitor MPP + was applied to 
FBX07 knockdown or over-expressing cells for 24 hours. MTT 
assay was performed to evaluate the cell viability. As shown in 
Fig. 6A, the knock down of FBX07 was confirmed as 76% 
(/> = 0.005-<0.001) and 81% (p = 0.004-0.001) of FBX07 were 
seen in FBX07 siRNA-transfected HEK-293T and SH-SY5Y 
cells, respectively, as compared to control siRNA-transfected 
(106%~ 1 10%) or untransfected (100%) cells. Fig. 6B examined if 
suppression or over-expression of FBX07 affects cell survival 
upon MPP + treatment. With 300 \lM MPP + treatment, the cell 
viability of HEK-293T cells expressing FBX07 siRNA was 61% 



of the untreated cells, which is significandy lower than that of 
untransfected (69%, p = 0.001) or control siRNA transfected (68%, 
p — 0.033) HEK-293T cells. A similar result was obtained in SH- 
SY5Y cells treated with 2 mM MPP + , as suppression of FBX07 
expression also significantly decreased the cell viability (64%) as 
compared to untransfected (73%, p = 0.024) or control siRNA 
transfected (74%, p= 0.018) cells. Conversely, significantly in- 
creased cell viability was observed with MPP + treatment in 
FBX07-EGFP over-expressing cells compared to untransfected 
cells (HEK-293T: 73% vs. 69%, ^ = 0.013; SH-SY5Y: 81% vs. 
73%, p = 0.037). 

SH-SY5Y cell model 

To test the effect of Cys52 on neuronal phenotype, we 
constructed Flp-In SH-SY5Y cells with Tyr52 or Cys52 
FBX07-EGFP expression in an inducible fashion. Immunoblot 
analysis showed that the FBX07 protein level was significandy 
increased in Cys52 cells as compared to that of Tyr52 cells after 
induction with doxycycline (+ Dox) for 2 days (128%, p = 0.042) 
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Table 1. Primers and conditions for PCR amplification 


of FBX07 cDNA and genomic DNA. 






Product/ 




Anneal (°C)/ 


RFLP enzyme 


Test (amplified region) 


MgCI 2 (mM) 


(fragment, bp) 


cDNA sequencing 


F: CTCTTTCCCCGTTTCGCC 


58/1.5 


1955 


R: GGAGAACCAAGAGCAGGGAGA 


pEGFP-N1-FBX07 cDNA cloning 


F: AAGCTTCTCTTTCCCCGTTTCGCCTCAG 


62/1 .5 


1727 


R: ACCGGTGGCATGAATGACAGCCGGCC 


pcDNA3.1/V5-His-FBX07 cloning 


F: AAGCTTCTCTTTCCCCGTTTCGCCTCAG 


62/1 .5 


1726 


R: CTCGAGACATGCATGACAGCCGGCCATC 


Y52C (TAC/TGC) 


F: AGGCTGAGGCAGGAGGAT TG 


58/1.5 


Psfl: CTGCAG 


R: CTCCAGTGAGGGGATCCCTG 




(240/222, 18) 


Ml 151 (ATG/ATA) 


F: TCACTGGAGATGAAGAGACC 


58/1.5 


Pael: GCATGC 


R: AATCGCTTGAACTGGGAGGC 




(300, 179/479) 


The Pst\ restriction site was created by PCR using a mismatch primer. For Y52C and Ml 1 51 amplification, the underlines 


in the primer sequence and enzyme recognition 


site indicate the mismatch nucleotide and polymorphic site, respectively. For cDNA cloning, the underlines in the primer sequence indicate the introduced H/ndlll, Age\ 



and Xho\ restriction sites. 
doi:1 0.1 371 /journal.pone.01 01 392.t001 



(Fig. 7A). Compared to the non-induced cells (- Dox), TRAF2 
protein level was significantly decreased in both Tyr52 (100% vs. 
120%,/;= 0.024) and Cys52 (87% vs. 128%, p = 0.022) FBX07- 
EGFP expressed cells. The difference in TRAF2 abundance 
between Tyr52 and Cys52 FBX07-EGFP expressed cells was also 
significant (100% vs. 87%, ^ = 0.042). These FBX07 ceUs were 
induced for differentiation with retinoic acid [16] for 7 to 21 days. 
Representative fluorescence microscopy images of cells differen- 
tiated for 21 days are shown in Fig. 7B. Significantly more total 
outgrowth in Cys52 cells was observed compared to Tyr52 cells 
after differentiation for 7-21 days (131-165%, p = 0.014-<0.001). 

Discussion 

Up to now, only four different types of FBX07 mutations 
(T22M, R378G, R498X, IVS7+1G/T) have been reported to be 
responsible for parkinsonian-pyramidal disease, which has been 
designated as the cause of PARK 15 [4,5,6]. Recently, two 
missense substitutions (p.Ile87Thr and p.Asp328Arg) were found 
in two EOPD patients in Taiwan [7]. We did not detect any 
mutation in our EOPD patients, which is compatible with the 
previous results of studies showing rare pathogenic mutations of 
the FBX07 gene in typical PD patients of Chinese ethnicity [1,7]. 
Nevertheless, our method using cDNA sequencing may miss some 
mutations in the non-coding and regulatory elements of the 
FBX07 gene. Also variants which result in a truncated protein or 
nonsense mediated decay would not be identified. 

Although the role of FBX07 in neurons is still not known, 
because it is a part of SCF ubiquitin ligase complex, its function in 
the ubiquitin-mediated protein degradation is implicated [9]. 
Impaired ubiquitin-mediated protein degradation has been found 
in sporadic PD [17], autosomal dominant PARK1 [18], and 
autosomal recessive PARK2 [19]. Therefore, it is postulated that 
FBX07 mutations may compromise the ubiquitin-proteasome 



function and cause neuronal dysfunction in PD. In this study, we 
showed that the Y52C AG genotype or G allele of FBX07 
conferred a reduced susceptibility to Chinese PD when the data 
from our study and Luo's [1] were combined. However, as Y52C 
G allele is rare and the difference in frequency is small (0.4% in 
PD and 1.2% in controls), the genetic evidence is very limited and 
the findings may be due to chance alone. While FBX07 mutations 
typically causes autosomal recessive parkinsonism with pyramidal 
tract signs, the Y52C AG genotype or G allele in this study 
appears to provide a protective effect in a dominant mode. T22M, 
R378G, R498X mutations resulted in decreased stability of 
FBX07 protein and T22M caused loss of its nuclear activity, 
both of which may jeopardize the neuronal function [15]. Using 
the program SWISS-MODEL, we showed that the Cys52 formed 
H-bond with the Asp54 causing decreased local energy, which 
may increase the stability of the protein. The increased stability of 
FBX07 protein consequent to Cys52 was further confirmed by 
cycloheximide chase experiment (Fig. 2). These results suggest that 
Cys52 may play a protective role in PD via increasing stability of 
FBX07 protein, which is in contrast to the decreased stability and 
loss of function caused by T22M, R378G, and R498X [15]. 
Indeed, knock down of FBX07 leading to dopaminergic neuronal 
dysfunction and cell loss and locomotor deficits that are improved 
by apomorphine, have been shown in zebrafish [20] . The study in 
zebrafish further suggests that FBX07 plays an important role in 
the development of dopaminergic neurons and its loss of function 
caused by mutations may be responsible for the phenotype of PD. 

FBX07 can interact with three proteins including hepatoma 
upregulated protein, cIAPl, and the proteasome inhibitor protein 
PI31 [1 1,21,22]. FBX07 was also reported to enhance activity of 
cyclin D/ cdk6 that plays an important role in regulating neuronal 
death processes [23]. How these interactions may contribute to the 
pathogenesis of PD remains to be clarified. 
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Table 2. Genotype and allele distributions and association analysis. 







No. (%) 




Odds ratio (95% CI) 


p-Value 




PD 


Control 






Genotype/allele 


Y52C (TAC/TGC) 


AA 


512 (99.2) 


504 (97.7) 


1.00 




AG 


4 (0.8) 


12 (2.3) 


0.33 (0.09-0.95) 


0.055 


Locus total p 




0.044 






A 


1028 (99.6) 


1020 (98.8) 


1.00 




G 


4 (0.4) 


12 (1.2) 


0.33 (0.09-0.95) 


0.056 


Locus total p 




0.045 






Ml 151 (ATG/ATA) 


AA 


258 (50.0) 


258 (50.0) 


1.00 




AG 


220 (42.6) 


221 (42.8) 


1.00 (0.77-1.28) 


0.972 


GG 


38 (7.4) 


37 (7.2) 


1.03 (0.63-1.67) 


0.914 


Locus total p 




0.992 






A 


736 (71.3) 


737 (71.4) 


1.00 




G 


296 (28.7) 


295 (28.6) 


1.00 (0.83-1.22) 


0.961 


Locus total p 




0.961 






Y52C (TAC/TGC) (Luo et al., 2010) 


AA 


1 33 (98.5) 


192 (96.0) 


1.00 




AG 


2 (1.5) 


8 (4.0) 


0.36 (0.05-1.47) 


0.202 


Locus total p 




0.184 






A 


268 (99.3) 


392 (98.0) 


1.00 




G 


2 (0.7) 


8 (2.0) 


0.37 (0.06-1.47) 


0.201 


Locus total p 




0.187 






Y52C Combined [Taiwan + China (Luo et al, 2010)] 


AA 


645 (99.1) 


696 (97.2) 


1.00 




AG 


6 (0.9) 


20 (2.8) 


0.32 (0.12-077) 


0.016 


Locus total p 




0.011 






A 


1 296 (99.5) 


1412 (98.6) 


1.00 




G 


6 (0.5) 


20 (1.4) 


0.33 (0.12-0.77) 


0.017 


Locus total p 




0.012 







Odds ratios were calculated by comparing each value with the major common genotype or allele. 
doi:1 0.1 371 /joumal.pone.01 01 392.t002 



More emerging evidence has suggested that neuroinflammation 
is involved in the pathogenesis of PD through inflammatory 
mediators such as TNFot, nitric oxide (NO), IL-6, and IL- 1 (3 [13]. 
NF-kB activation is required for all of these inflammatory 
mediators to be produced by microglial cells [14]. Recently, few 
NF-kB inhibitors have been applied to the therapeutic approaches 
of several chronic inflammatory diseases including PD [13]. 
Interestingly, FBX07 protein has been shown to be a negative 
regulator of NF-kB signaling pathway, through binding to and 
ubiquitinating TRAF2 and cIAPl, which would lead to decreased 
RIP1 ubiquitination and NF-kB activity [12]. To further 
investigate the role of FBX07 in NF-kB signaling, Cys52 variant 
was expressed in HEK-293T and SH-SY5Y cells and TRAF2 
(NF-kB signaling protein) was examined in this study. As shown in 
Fig. 4-6, structurally stable Cys52 FBX07 facilitated the 
degradation of TRAF2 protein via increasing TRAF2 ubiquina- 
tion, which may theoretically further lower the NF-kB activity. 
Although the NF-kB activity was not evaluated in this study, the 
results provide some evidence that Cys52 variant may play a 



protective role in PD pathogenesis. Whether other NF-kB 
signaling proteins such as RIP1 are also regulated by FBX07 
needs to be investigated in future studies. 

In conclusion, gene mutations may be rare in Chinese early- 
onset Parkinsonism patients. We have shown that Y52C 
polymorphism of FBX07 may contribute to reduced PD 
susceptibility in Chinese population. However, additional case- 
control studies are needed to establish whether FBX07 variants 
truly play a role in PD. 

Materials and Methods 

Ethics statement 

This study was performed according to a protocol approved by 
the Institutional Review Board of Chang Gung Memorial 
Hospital, and all examinations were performed after obtaining 
written informed consents. 
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Figure 2. Expression of FBX07 protein in HEK-293T cells. (A) Confocal microscopy examination of Tyr52 and Cys52 FBX07-EGFP protein in 
HEK-293T cells. HEK-293T cells were transfected with Tyr52 and Cys52 FBX07-EGFP for two days and stained with DAPI to detect nuclei. The green 
(FBX07 fusion protein) to blue (nuclei staining) fluorescence signal was quantified and shown below. (B) Western blot analysis of pEGFP-N1 vector 
(Vec) and Tyr52 (WT) and Cys52 FBX07-EGFP transfected cells using FBX07 antibody and anti-tubulin and anti-neomycin antibodies as loading and 
transfection controls, respectively. (C) Tyr52 and Cys52 FBX07-EGFP proteins in transfected cells were analyzed after blocking of new protein 
synthesis in a cycloheximide (200 |ig/ml) chase experiment (cells were harvested after 6, 12, 24, 36, and 48 hr). Data are represented as the means ± 
S.D. of three separate experiments. 
doi:1 0.1 371 /journal.pone.01 01 392.g002 



Subjects 

A total of 516 unrelated Taiwanese PD subjects (45.0% females) 
were recruited from the neurology clinics of Chang Gung 
Memorial Hospital (CGMH). All patients were diagnosed by 
two neurologists specialized in movement disorders (Y.-R. Wu and 
C.-M. Chen) with probable idiopathic PD according to the 
published criteria [24], which includes substantial and sustained 
response to levodopa or a dopamine agonist. Subjects with prior 
history of multiple cerebrovascular events or other causes of 
parkinsonian symptoms (e.g. brain injury or tumor, encephalitis, 
antipsychotic medication) were excluded. The mean age at onset 
(AAO) of PD was 62.0± 11.5 years, ranging between 19 and 93 
years. For juvenile PD patients (AAO<50), mutations in the 
Parkin, PLNK1, DJ-1, ATP13A2 and G2019S LRRK2 were excluded 
([25,26,27] and unpublished results). A group of 516 normal 
controls without neurodegenerative diseases were recruited from 
the same ethnic community. Control subjects (50.2% females) had 
mean age at examination of 60. 9± 12.3 years, ranging between 20 



and 92 years. All examinations were performed after obtaining 
written informed consent from patients and control individuals. 
This study was approved by the Institutional Review Board of 
CGMH. 

Genetic analysis 

Genomic DNA was extracted from peripheral blood leucocytes 
using the standard protocols. For PD patients with onset ^50 
(n = 80, mean age at onset 43. 7 ±0.7 years, 33.7% females), RNA 
was extracted using PAXgene Blood RNA Kit (PreAnalytiX). The 
RNA was DNase (Stratagene) treated, quantified, and reverse- 
transcribed to cDNA using High Capacity cDNA Reverse 
Transcription Kit (Applied Biosystems). Using polymerase chain 
reaction (PCR) with designed primers and conditions (Table 1), 
the 1955-bp amplified FBX07 cDNA was gel purified and 
sequenced directly using the ABI PRISM 3 1 30 Genetic Analyzer 
(Applied Biosystems). The identified Y52C and Ml 151 variants 
were verified by genomic DNA PCR and sequencing. For 
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Figure 3. The homology models of Tyr52 (WT) FBX07 and Cys52 (Y52C). Protein models were shown in secondary representation. Tyr52 
and Cys52 residue were shown in amino acid structures. The hydrogen bond interaction was presented in yellow dash line. 
doi:1 0.1 371 /journal.pone.01 01 392.g003 



population screening, the Y52C and Ml 151 were examined using 
the Pstl and Pael (gain of sites) restriction enzymes, respectively 
(Table 1). The digested PCR products were visualized with 
ethidium bromide after electrophoresis in 2.2% or 1.6% agarose 
gel. 



FBX07 cDNA constructs 

Using the designed primers to remove translation termination 
codon (Table 1), the full-length FBX07 cDNA fragments from an 
individual heterozygous for Cys52 were cloned into pGEM-T 
Easy vector (Promega) and sequenced. The 1.7 kb Hindlll (added 
in the forward primer)-4gd (added in the reverse primer) 
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Figure 4. Western blot analysis of NF-kB signaling pathway protein TRAF2. HEK-293T cells were tranfected with Tyr52 or Cys52 FBX07 
construct and harvested after two days. Transfected FBX07-EGFP or endogenous TRAF2 protein was detected using specific antibodies and anti-actin 
and anti-neomycin antibodies as the loading and the transfection control, respectively. Data are represented as the means ± S.D. of three separate 
experiments. 

doi:1 0.1 371 /journal.pone.01 01 392.g004 
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Figure 5. Co-immunoprecipitation of Tyr52 or Cys52 FBX07-EGFP protein and ubiquitination of TRAF2. (A) HEK-293T cells were 
transiently transfected with Tyr52 or Cys52 FBX07-EGFP construct. After 48 h, cell lysates were prepared (Input, left panel) and immunoprecipitations 
(IP, right panel) were performed with anti-TRAF2 antibody. Normal IgG was used as a negative control for IP. Cell lysates and immunoprecipitates 
were analyzed with anti-EGFP, anti-TRAF2, anti-ubiquitin or anti-actin antibody. (B) Quantification of immunoprecipitated FBX07-EGFP, TRAF2 and 
ubiquitin in HEK-293T cells transiently transfected with Tyr52 or Cys52 FBX07-EGFP construct for 2 days. Data are represented as the means ± S.D. of 
three separate experiments. 
doi:1 0.1 371 /journal.pone.01 01 392.g005 



fragments were removed from pGEM-T Easy vector and ligated 
into the corresponding sites of pEGFP-Nl (Clontech) to generate 
Tyr52 and Cys52 FBX07 cDNA in-frame fused to the EGFP gene. 
The resulting EGFP-tagged FBX07 constructs were used in 
transient expression studies for confocal microscopy examination, 
FBX07 stability and anti-TRAF2 co-immunoprecipitation. Ad- 
ditionally, the Hindlll-Xhol fragments containing FBX07 were 
ligated into pcDNA3.1/V5-His (Invitrogen) to generate Tyr52 and 
Cys52 FBX07 cDNA in-frame fused to the V5-His for Western 
blot analysis of NF-kB signaling pathway protein TRAF2. 

Cell cultivation and transfection 

Human embryonic kidney (HEK)-293T (ATCC No. CRL- 
11268) cells were cultivated in Dulbecco's modified Eagle's 
medium containing 10% fetal bovine serum (FBS) in a 37°C 
humidified incubator with a 5% CO2 atmosphere. Cells were 
plated into 6-well (6xl0 5 /well) dishes, grown for 20 hr and 
transfected by the lipofection method (GibcoBRL) with EGFP- 
tagged FBX07 constructs (4 u,g/well). The cells were grown for 



48 hr for the protein studies. To evaluate the stability of FBX07 
protein, protein synthesis inhibitor cycloheximide (200 u,g/ ml) was 
added 24 hr after transfection for 0, 6, 12, 24, 36, and 48 hr 
before protein preparation. For immunoprecipitation studies, 
proteasome inhibitor MG-132 (5 uM) was added 24 hr after 
transfection for 24 hr before protein preparation. 

Confocal microscopy examination 

For visualizing intracellular FBX07-EGFP protein, transfected 
cells on coverslips were stained with 4'-6-diamidino-2-phenylin- 
dole (DAPI) to detect nuclei. The stained cells were examined for 
dual fluorescent imaging using a Leica TCS confocal laser 
scanning microscope. 

siRNA transfection and cell viability assay 

Human neuroblastoma SH-SY5Y cells (ATCC No. CRL-2266) 
were maintained in DMEM F12 supplemented with 10% FBS at 
37°C in an atmosphere containing 5% C0 2 . For FBX07 
knockdown, siRNA specifically targeting FBXQ7 (J-0 13606-06, 
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Figure 6. Toxic effects of MPP + in FBX07 suppressed or over-expressed cells. (A) HEK-293T and SH-SY5Y cells were transfected with FBX07- 
specific or control siRNA. After 48 h, cell lysates were prepared and Western blot analysis was performed using FBX07 and anti-actin (as loading 
control) antibodies. (B) HEK-293T and SH-SY5Y cells were transfected with siRNA (FBX07-specific or control), FBX07 cDNA or not (-). After 24 hr, cells 
were treated with MPP + (300 nM for HEK-293T or 2 mM for SH-SY5Y) and toxic effects of MPP + were monitored at 24 hr by MTT assay. 
doi:1 0.1 371 /journal.pone.01 01 392.g006 



Thermo Scientific) was transfected into in HEK-293T (by 
lipofection) or SH-SY5Y cells (by 4D-Nucleofector System, 
Lonza). For transfection of SH-SY5Y, cells (2.5 x 10 r 748-well) 
were centrifuged at 90 xg for 10 min at room temperature. The 
cell pellets were resuspended in 4D-Nucleofector Solution 
(100 U.L) with siRNA (100 nM) or FBX07-EGFP (250 ng) and 
CA-137 appropriate program was selected. The transfected cells 
were plated in 48-well plate and grown for 48 hr for protein 
analysis. Alternatively, after transfection for one day, cells were 
treated with freshly prepared mitochondrial inhibitor MPP + (1- 
methyl-4-phenylpyridinium ion, 300 |J,M for HEK-293T or 2 mM 
for SH-SY5Y cells). After one day, 20 uL MTT [3-(4,5- 
dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide, 5 mg/mL 
in PBS, Sigma] was added to cells and incubated for 2 hr. The 
absorbance of the insoluble purple formazan product was 
measured at 570 nm by a Bio-Tek u.Quant Universal Microplate 
Spectrophotometer. 



Western blot analysis 

Cells were lysed in hypotonic buffer (20 mM HEPES pH 7.4, 
1 mM MgCl 2 , 10 mM KC1, 1 mM DTT, 1 mM EDTA pH 8.0) 
containing the protease inhibitor mixture (Sigma). After sonication 
and sitting on ice for 20 min, the lysates were centrifuged at 
14,000 xg for 30 min at 4°C. Protein concentrations were 
determined using the Bio-Rad protein assay kit, with albumin as 
standards. Total proteins (25 u,g) were electrophoresed on 10% 
SDS-polyacrylamide gel and transferred onto nitrocellulose 
membrane (Schleicher and Schuell) by reverse electrophoresis. 
After being blocked, the membrane was stained with anti-FBX07 
(1:3000 dilution, Abnova), anti-TRAF2 (1:500 dilution, Santa 
Cruz), anti-neomycin (1:1000 dilution, Millipore), anti-tubulin 
(1:10000 dilution, GeneTex), anti-GAPDH (1:1000 dilution, 
MDBio), anti-actin (1:10000 dilution, Millipore), anti-EGFP 
(1:200 dilution, Santa Cruz), or anti-ubiquitin (1:1000 dilution, 
Dako) antibody. The immune complexes were detected using 
horseradish peroxidase-conjugated goat anti-mouse (Jackson 
ImmunoResearch) or goat anti-rabbit (Rochland) IgG antibody 
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Figure 7. Flp-ln SH-SY5Y cells with induced FBX07 expression and neuronal phenotype. (A) Western blot analysis of FBX07-EGFP protein 
level using FBX07 antibody and anti-actin antibody as loading control after 2 days induction with doxycycline (+ Dox) or not (- Dox). (B) 
Representative microscopic images of neuronal differentiated Tyr52 and Cys52 cells (for 21 days). Nuclei were counterstained with DAPI (blue). 
Neuronal total outgrowth was quantified in Tyr52 and Cys52 cells with induced differentiation for 7-21 days. Data are represented as the means ± 
S.D. of three separate experiments. 
doi:1 0.1 371 /journal.pone.01 01 392.g007 



(1:10000 dilution) and Immobilon Western Chemiluminescent 
HRP substrate (Millipore). 

Immunoprecipitation 

Total protein from MG-132-treated FBX07-EGFP-transfected 
cells was prepared using IP lysis buffer (25 mM Tris-HCl pH 7.4, 
150 mM NaCl, 1% NP-40, 1 mM EDTA, 5% glycerol) contain- 
ing the protease inhibitor mixture. After quantification, proteins 
(100 |xg) were immunoprecipitated with mouse TRAF2 primary 
antibody (2 Jig, Santa Cruz) conjugated to protein G beads 
(Millipore). The beads-proteins-antibody mixtures were washed 
three times with IP lysis buffer and the immunoprecipitated 
proteins were eluted by 1 x SDS sample loading buffer (50 mM 
Tris pH 6.8, 2% SDS, 10% glycerol, 2.5% P-mercaptoethanol, 
0.005% bromophenolblue), separated on 10% SDS-polyacryl- 
amide gel, transferred onto nitrocellulose membrane and probed 
with indicated antibodies as described. 

FBX07 SH-SY5Y cell lines generation 

The SH-SY5Y-derived Flp-ln host cells [28] and Flp-ln T-REx 
System (Invitrogen) was used to generate stably induced SH-SY5Y 
cell lines exhibiting tetracycline-inducible expression of Tyr52 and 
Cys52 FBX07. Briefly, the SH-SY5Y host ceUs were co- 
transfected with pOG44 plasmid (constitutively expressed the Flp 
recombinase) and pcDNA5/FRT/TO-FBX07-EGFP plasmid 
according to the supplier's instructions. These cell lines were 
grown in medium containing 5 jxg/ml blasticidin and 100 u.g/ml 
hygromycin. Doxycycline (dox, 5 |Xg/ml) was added to induce 
EGFP-tagged FBX07 expression for two days. The proteins were 
prepared for Western blotting using antibody to FBX07 or actin 



as described. Neuronal phenotypes were examined after induced 
differentiation with retinoid acid (10 uM) and induced expression 
of FBX07 for 7 to 21 days. The morphologic differentiation of 
Tyr52 and Cys52 SH-SY5Y cells including total outgrowth, 
processes, and branches was assessed by using Metamorph 
microscopy automation and image analysis software (Molecular 
Devices). 

Statistical analysis 

The genotype frequency data and the expected genotypic 
frequency under random mating were computed and Chi-square 
tested for Hardy-Weinberg equilibrium using standardized 
formula. The genotype and allele association analysis was carried 
out using the Chi-square test. Odds ratios with 95% confidence 
intervals (95% CI) were calculated to test association between 
genotype/allele and disease. Given the observed Y52C G allele 
frequency of 0.0078 (0.0095, in combined data from Taiwan and 
China) and a total of 1032 (1367, combined subjects from Taiwan 
and China) subjects in the present study, at significance level of 
0.05, we had power greater than 0.8 to identify an association 
when the allele genetic effect size was greater than 3. 1 (2.5). Given 
the observed Ml 151 G allele frequency of 0.286 and 1032 subjects 
in the present study, at significance level of 0.05, we had power 
greater than 0.8 to identify an association when the allele genetic 
effect was greater than 1.5. 

For statistical analysis of microscopy images, immunoblots and 
cell viability assays, data were expressed as the means ± standard 
deviation (SD). Three independent experiments were performed 
and non-categorical variables were compared using the Student's 
/-test. All /(-values were two-tailed, with values of p<0.05 being 
considered significant. 
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Homology modeling 

We modeled the three dimensional structures of the Tyr52 and 
Cys52 FBX07 proteins by comparative methods and energy 
minimization using the program SWISS-MODEL [29]. The 2.9-A 
coordinate set for the crystal structure of human UBC protein 
(PDB code 2ZVO, chain A) served as the template for modeling 
the residue 1-79 of human FBX07. The energy computation was 
done with the GROMOS96 [30] implementation of Swiss- 
Pdb Viewer. The resulting FBX07 three-dimensional models were 
manipulated and rendered in PyMOL (The PyMOL Molecular 
Graphics System, Version 1.2r3pre, Schrodinger, LLC). 
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